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Abstract. We calculate the perturbed action, at second and third order, for a massive 
three-form field minimally coupled to gravity, and use it to explore the observational predic- 
tions of three-form inflation. One intriguing result is that the value of the spectral index is 
nearly independent of the three-form potential, being fixed solely by the number of e-folds 
of inflation, with n s = 0.97 for the canonical number of 60. Considering the bispectrum, we 
employ standard techniques to give explicit results for two models, one of which produces a 
large non-Gaussianity. Finally, we confirm our results by employing a duality relating the 
three-form theory to a non-canonical scalar field theory and explicitly re-computing results 
in this dual picture. 
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1 Background 



Cosmological inflation [1—6], typically assumed to be driven by one or more scalar fields, 
offers a compelling explanation for the origin of structure in our universe (see [7-10] for 
reviews). But if inflation is indeed the source of cosmic structure, the field or fields driving 
this inflationary period have so far not been identified, and a great many candidates have 
been suggested, including scalar fields with non-canonical kinetic terms [11—13]. 

During inflation, quantum fluctuations are promoted to classical perturbations, and the 
statistics of these perturbations are probed by observations (for example [14, 15]). In recent 
years, interest in the statistics beyond the two point function, which is parametrised by the 
power spectrum, has intensified. In particular, numerous investigations of the three-point 
function, parametrised by the bispectrum, have been undertaken [16-18]. The bispectrum is 
potentially a powerful diagnostic, for example it may allow non-canonical single field infla- 
tion, which is capable of producing an observable bispectrum [18], to be distinguished from 
canonical inflation, for which the bispectrum is typically too small ever to be observed [16]. 

1.1 Three- form inflation 

In addition to scalar fields, other ways of driving inflation have been proposed. These include 
inflation driven by vector fields (for example [19-21]) and higher form fields [22-27]. For such 
models, the techniques needed to make observational predictions are less well explored than 
for scalar field models. In this paper, our interest lies in the recently proposed inflationary 
scenarios involving a massive three-form field. Induced potentials for three-forms were studied 
in the context of string theory in Refs.[28, 29]. Inflationary scenarios with three- forms were 
first investigated by Germani & Kehagias [22, 25], who focused on a non-minimally coupled 
three-form, with the couplings chosen to make the theory exactly equivalent to a canonical 
and minimally coupled scalar field theory. The behaviour of gravitational waves in this 
models was studied in Ref.[30]. Later Koivisto & Nunes [23, 24] considered a minimally 
coupled massive three-form, which has a novel and much richer associated phenomenology. 
They found that if particular conditions on the potential were satisfied, the model was ghost- 
free, and moreover could support inflation. More recently, Koivisto & Urban [31] studied the 
non-Gaussian signatures resulting from the magnetic fields generated from the coupling of a 
three-form to electromagnetism. 

Here we follow Koivisto & Nunes and consider the action for a minimally coupled, 
canonical three-form 



where 8nG = 1, F = VA is a four- form, A a three- form, and we have adopted a compact 
notation in which indices are suppressed once the valance of an object is specified, squaring a 
quantity denotes contraction of all indices (for example F 2 = F pupu F^ upa = V p A upa V fl A up ' T ), 
and a dot (•) denotes contraction on the first index (for example V • A = V^Apa). Here and 
throughout Greek indices label space-time dimensions and Roman indices spatial dimensions. 

In the context of inflationary cosmology, the unperturbed three-form must be consistent 
with the universal symmetries of homogeneity and isotropy. The background three-form can 
then be written in the form 




(1.1) 



Aijk = a(t) 3 e ijk x(t) 



(1.2) 
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where a(t) is the scale factor of the universe as a function of cosmic time, and all other 
elements of A are zero. Hence the unperturbed three- form satisfies A 2 = 6x 2 • The dynamics 
of the universe is then governed by the behaviour of the scalar quantity x(t) which is directly 
related to the three- form. As emphasised elsewhere [23, 24], when written in terms of this 
scalar quantity the equations of motion which govern the behaviour of the universe, and the 
role of the three-form potential, are straightforward to interpret. In particular, adopting 
units where 8irG = 1, the Friedmann and Raychaudhuri equations are given by 



respectively. 

On the other hand, as we will review in §3, a series of dualities exist which relate the 
action (1.1) to the action for a massive vector field, and in turn to the action for a massive, 
and in general non-canonical, scalar field of the form 



where 4> is a scalar field and X = — V^V^. In principle, therefore, if one wished to 
understand the dynamics and observational predictions of a particular three-form inflationary 
model, one could move to working with the dual non-canonical scalar field model, for which 
suitable techniques are already readily available [12, 16-18]. Unfortunately the mapping 
between theories is rather complicated, and the analytical form of the dual non-canonical 
scalar field theory is often rather unwieldy or may not even exist, even for fairly simple 
three-form potentials. 

With this in mind, Koivisto & Nunes [23, 24] advocated working directly with the 
three-form. Utilising the Einstein field equations, they derived the background equations 
of motion, and perturbing both the three-form and the metric, the first order perturbation 
equations. Finally, they formed the Sasaki-Mukhanov equation for the curvature perturbation 
on comoving hypersurfaces, £. This equation follows from the second order perturbed action 
written purely in terms of £, which they thus implied. In addition they were able to address 
questions such as the stability of the theory, and the power spectrum produced by three- form 
inflationary models. Later De Felice, Karwan & Wongjun [27] directly perturbed the action, 
and, following the standard approach, rewrote the result purely in terms of the curvature 
perturbation. In this way they obtained similar results, though curiously their perturbed 
action did not match exactly that of Koivisto & Nunes. 

1.2 Overview 

In this paper, we initially follow Koivisto & Nunes by working with the three-form theory 
itself rather than any dual. In contrast to their study, however, we manipulate the action 
directly, perturbing it to second order and writing the result purely in terms of £. Doing 
so we find agreement with their result. Then we extend existing work by deriving the third 
order action, which is required to calculate the bispectrum of the curvature perturbation 
produced during inflation. We also explain why the two previously derived second order 

1 Note that in Ref. [24], the three form quantity was denoted by X not Xi but here we reserve X for use 
when we discuss non-canonical scalar field models to match notation commonly used in the literature. 



3H 2 = - (x + 3H X ) 2 + V( X 2 ) 
H = -\v x x, 



(1.3) 




(1.4) 
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actions, by Koivisto & Nunes and De Felice, Karwan & Wongjun [27] respectively, differed 
slightly. These calculations are covered in §2. In §3, we review the dualities which relate 
the three-form theory to a non-canonical scalar-field theory. We then show that despite 
the complexity of the mapping between the theories, we may use this duality to provide an 
alternative derivation of our second and third order action, providing confirmation of our 
results. Finally, we employ our results to calculate the bispectrum in two explicit examples 
of three-form inflation, one of which produces a large non-Gaussianity, in §4. We conclude 
in §5. 

2 Perturbations and observables 

In order to make observationally testable predictions for inflation driven by a three-form field, 
we must understand the statistics of perturbations produced during the inflationary era. In 
particular, it is usual to focus on £ [32, 33], the co-moving curvature perturbation, which is 
conserved on super-horizon scales if the dynamics is adiabatic [34, 35]. This property helps 
to connect the statistics of the Fourier modes of £ soon after they cross the horizon with 
observations of later epochs. 

To this end, we wish to perturb the action for the three-form theory (see [36-38] for 
reviews of cosmological perturbation theory). The aim is to determine the form of the 
autonomous action for £ in the setting of three form-inflation. Then by following the standard 
procedure of quantising £ and fixing initial conditions far inside the cosmological horizon, the 
same-time correlation functions of £, produced during inflation, can be calculated using the 
In-In formalism. The second order action is required to calculate the two-point correlation 
function, parametrised by the power spectrum (see for example [7] for a review), while the 
three-point correlation function, parametrised by the bispectrum, requires knowledge of the 
third order interactions in the theory, and hence the third order action [16-18]. 

The procedure to calculate the action is to perturb the metric and the three-form, and 
insert the result into Eq. (1.1). Using various constraint equations, it is then possible to 
write the result purely in terms of only the perturbed quantity £, and background quantities. 
The action can then be separated out order by order. We now present the details of this 
calculation. 

2.1 The perturbed action for three form inflation 

Following the approach of Maldacena [16], we consider the Arnowitt-Deser-Misner form of 
the metric 

ds 2 = -iV 2 di 2 + hij (dx i + iVMt) (dx j + N'dt) , (2.1) 

where ./V is the lapse function, Ni the shift vector, and the spatial metric. In this notation 
the action for the three-form takes the form 

S =\J dtdx 3 VhN (r® - i_F 2 - 2V{A 2 )J + dtdx^VhN' 1 (EijE ij - F 2 ) , (2.2) 

where h = dethij, is the Ricci scalar associated with the spatial metric, and F^ = 
1/2/ijj — Nuu\ is the extrinsic curvature. Considering only scalar perturbations, when per- 
turbed about an FRW metric the spatial metric takes the form 

hij = a 2 e 2 % , (2.3) 
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where we have used the spatial gauge freedom available to make the spatial metric diagonal, 
and where £ denotes the curvature perturbation. The lapse function and shift vector are also 
perturbed to give 

Ni = il> ti + Ni, N = 1 + 5, (2.4) 

where N% i = and if) and a can be written as an expansion in powers of £ as d = a.\ + a.<i 
and ip = ipi + -02 • 

We must also account for perturbations in the matter sourcing the Universe's evolution, 
which in this case is the three-form field. The most general form of the perturbed 3-form is 

A ij = a(t)e ijk (a k + a k ), A ijk = a(t) 3 e ijk (X(t) + a ) . (2.5) 

The vector perturbation a k together with the vector perturbation in the metric (which we 
did not write down) decouple from scalar perturbations in the second order action and also 
decay as usual [24]. At higher order they can also be neglected since our interest here is in 
the interaction terms which contribute to the scalar-scalar-scalar three-point function. The 
metric tensor perturbation also decouples in the second order action and can be neglected at 
higher order for the same reason. There are two scalar degrees of freedom, a and ao, which 
just like the metric perturbations can be expanded in powers of £. 

All the scalar perturbed quantities for the both the metric and the matter are then 
substituted into Eq. (2.2), which is expanded to third order. 

2.2 Fixing the temporal gauge 

Thus far in our calculation we have only fixed the spatial part of the gauge freedom. At first 
order in perturbation theory, under a general gauge transformation, — > x u + where 
£^ = (£°>£i) the scalar three- form perturbations transform as 

a ^ ao - + IV 2 ? , a -> a - aX£ . (2.6) 

The freedom associated with £ is already fixed, but not that associated with £°, Considering 
the transformations above, it is clear that we could fix the gauge by setting ao to zero, which 
was the choice made by De Felice, Karwan & Wongjun [27] . This can also be done at second 
order. 

However, it is important to note that this does not correspond to the choice which is 
commonly made, for example by Maldacena [16], Seery & Lidsey [17] and Chen et al. [18], 
who select the comoving gauge, defined at linear order by the condition 

5T? = 0. (2.7) 

In models sourced by a scalar field, this is equivalent to setting 8(f) to zero, while in the 
present case it implies the following condition 

a = a 2 xV>- (2-8) 

Enforcing this condition at all orders fixes our temporal gauge freedom, and we can use this 
relation to replace a with if), or visa versa. We will later confirm explicitly that this choice 
is consistent with that made in work on scalar field inflationary models. 

Before moving on with our calculation, we note that in their study Koivisto & Nunes 
[23, 24] worked in the longitudinal gauge, but wrote down the action for the curvature 
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perturbation in the comoving gauge, £. The choice of De Felice, Karwan & Wongjun [27] 
selects a subtly different gauge, and explains why the second order actions in these two works 
do not coincide 2 . 

2.3 The second order action 

Moving forward with our calculation, we use the condition (2.8) to substitute for a in the 
action. We find a rather cumbersome expression containing terms in ip, ao, a and £. The first 
order part of the action is zero once the background equations (1.3) of motion are utilised 
(or equivalently varying the first order part of the action leads to the background equations 
of motion). Our aim now is to find substitutions for ip, ao and a in terms of C, and hence to 
obtain an action only in terms of £. As long as these substitutions follow from a constraint 
equation, and do not alter the order of the action, we do no harm by using them to write the 
action in terms of £ alone. First we focus on the second order part of the action. 

In the present case the simplest way to proceed is to utilise four constraint equations, 
which are not independent of each other, for the three quantities we wish to eliminate. The 
first is the equation of motion for N, or equivalently a 

3H(( -Ha)-~(x + 3H X )(a Q + 3Ha ) + \(x + 3#x) 2 (a + 3() 

-\v, x a + \xV x C + \a\{x + 3ff*)^ - ^(H^ + = 0. (2.9) 

The next is the equation of motion for N{ , or equivalently if) 

5 = |, (2.10) 

Then we can also use the equation which follows from the anti-symmetry of the four-form, 
and an equation which is most easily seen from the conservation of the energy momentum 
tensor. These are 

«o = 3 X C-^«i, (2.11) 
v ,xx 

d 2 

»o = a 2 X ^ + {x + 3H X )(a + 3() - 3Ha , (2.12) 

respectively. One can verify that all these equations are self-consistent both with each other 
and with the variation of the action with respect to a, if) and ao- 

We find that the substitutions which consistently satisfy the equations above are 

a=j ! , (2.13) 

ip = -jj+<J, (2.14) 
8 2 a = a 2 4c, (2.15) 



2 It is important to note, that there is no problem with picking a gauge other than the comoving gauge 
and then forming £ out of quantities in the chosen gauge at the end of the calculation. For example in Refs 
[16, 39] the consistency of perturbing the action in the fiat gauge and then calculating the statistics of £ was 
shown. 
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with e = —H/H 2 = xV,x/2H 2 and the speed of sound is given by [24] 

c 2 = (2.16) 

v ,x 

We note that these substitutions are only first order, while we might expect the second 
order action to include terms which come from 012 etc.. As highlighted by Maldacena [16], 
however, that these relations come from a constraint equation is advantageous, since it means 
that only the first order term in the substitutions for the perturbed quantities in terms of C, is 
required. The second order part of each term (for example a 2 ) will multiply the zeroth order 
constraint in the second order action. Looking ahead to the third order action, we note that 
even there we will only need the first order substitution. This is because the second order 
parts multiply the first order constraint in the third order action, and so disappear. While 
the third order part will multiply the zeroth order constraint in the third order action. That 
this occurs can be verified explicitly. 

The next step is to make all these substitutions in the second order part of the action 
so as to write it purely in terms of Q. After proceeding in this way, utilising the background 
equations (1.3) and integrating by parts in a suitable manner we arrive at the action 



S 2 = I dtd 6 x 

with 



a^C 2 -ae(d() 2 



(2.17) 



tt2 

S = =V- (2.18) 

We note that in this calculation the tp substitution is not actually used, as ip naturally cancels 
out of the second order action after the other substitutions. This second order action matches 
precisely that derived by Koivisto & Nunes [24]. 

2.4 The third order action 

Now we wish to isolate the third order part of the action in terms of £ alone. Once again we 
can make the substitutions above in the third order part. After numerous integrations by 
parts we arrive at 

S 3 = J Ltd 3 x[-eaad0 2 ~ a 3 (£ + 2A)-J + CC 2 

+ ( 3 ^ " jj ) [didj^did^ - 8^d 2 ij) - 2a- 1 d i ^diCd 2 ^ \ (2.19) 



where 



1 V%V XXX 
,xx 

and where for simplicity we have not explicitly substituted for i/j, but could do so using 
Eq. (2.14). 

At this juncture we pause to make an observation. The form of both the second and 
third order action, is precisely that found for a P(X, 4>) action of the form (1.4), except that 
in that case, c 2 and A are written in terms of derivatives of P(X, 4>) with respect to X [17]. 
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One might have expected this, given the formal equivalence of the theories. However, as we 
have noted previously, the mapping between them is rather complex, and working directly 
with the scalar dual for a three-form model with a given potential is often cumbersome or 
indeed impossible, if the analytic form of the inverse potential is not known 3 . We will see 
below, however, that using the duality in a more formal manner allows us to use the P(X, eft) 
description to recover Eqs. (2.16) and (2.20), effectively providing an alternative derivation 
for the perturbed three-form action to third order, and providing a powerful confirmation of 
our results thus far. We now proceed to that calculation. 

3 Dual theories 

In this section we review a number of dualities which inter-relate p-form theories in four 
space-time dimensions. These have been discussed elsewhere [22, 24-26], and our aim here 
is to employ them as a practical tool. To this end we show how one of these dualities can be 
used to confirm our earlier results for the perturbed action of a three-form theory, utilising 
previous results derived for a non-canonical scalar field theory by other authors [17, 18]. 

3.1 Equations of motion for the three-form 

We begin by considering a Palatini-type action, which treats the three-form A and the four- 
form F as independent variables, and which is equivalent to (1.1) up to boundary terms. The 
Lagrangian for the matter part of the action is given by 



Variations of the associated action relate A and F. The resulting equations of motion are 



where a prime denotes differentiation with respect to the argument in brackets. Integrating 
the middle term of C\ by parts inside the action, one finds 



which consequently shares the same equations of motion. The first equation of (3.2) now 
appears as a constraint equation and may be substituted back into £2 to confirm that we 
do indeed recover (1.1). We note that we are always free to perform such an integration by 
parts and this procedure leaves the dynamics of the theory invariant. Moreover, we are free 
to substitute constraint equations which do not change the order back into the action. 

3.2 Dual actions 

Our primary aim is to rewrite Lagrangians (3.1) and (3.3) in terms of the Hodge (*) dual 
fields to A and F. In appendix A, for clarity, we give some pedagogical detail regarding such 
duals. We recall that any p-form has a dual [d— p)-form, where d is the number of space-time 

3 This can be the case, for example, if the potential is a polynomial of degree five or above. 




(3.1) 



F = -4 [VA] 
V • F = -12 AY (A 2 ) 



(3.2) 




(3.3) 
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dimensions, four in our context. In particular, the three-form A and four-form F that make 
up (3.3) can be expressed in terms of their duals as 

(*F) = ^e a ^ s F a ^ 5 = $ F aM& = -e a ^ 

(*A) a = —^ea^sA^ 5 = B a Ap lS = -e a ^B a (3.4) 

The Hodge duals to F and A enable us to recast the original theory (1.1) into a scalar- 
vector description, with Lagrangian 

£ 3 = _l$ 2 _$V-S-y(-6 J B 2 ), (3.5) 

which follows from Lagrangian (3.3). The equations of motion for $ and can now be 
obtained either by varying (3.5) or equivalently by substituting Hodge duals into (3.2). They 
are 

= —V • B, 

V$ = -12BV'(-6B 2 ). (3.6) 

The first equation of motion now appears as a constraint equation with respect to £3 and can 
be substituted back to express our original theory (3.1) as a pure vector theory. Integrating 
the middle term of £3 by parts (or equivalently substituting A and F for their duals in 
£2), however, one finds that the converse is true, and the second equation of motion in 
(3.6) appears as a constraint and may be substituted into the action to remove the vector 
field in favour of the scalar (or equivalently its dual four form). With an eye on calculating 
inflationary observables, and in particular the 3-point correlation function, the scalar picture 
is particularly intriguing, since, as we have discussed, it opens up the possibility of using 
existing machinery for dealing with scalar field models of inflation. The final set of equivalent 
actions for a four-form, three-form, vector and scalar respectively are 

£iv(F,V -F) = -^-F 2 + 2A 2 (V ■ F)V' (A 2 (V ■ F)) - V '{A 2 (V • F)) , (3.7) 
48 

C IH (A,VA)) = -i [V^] 2 - V(A 2 ), (3.8) 

d(B,V-B) = ±(V-B) 2 -V(-6B 2 ), (3.9) 

V*) = -^$ 2 - 12B 2 (V<S>)V (-65 2 (V$)) - V (-65 2 (V$)) . (3.10) 

Here V (— 6-B 2 (V < ]?)) and V(A 2 (V • F)) indicate that the second equation of motion in (3.6) 
has to be used in order to express A 2 in terms of V • F and B 2 in terms of V<3?. Note that 
we cannot simultaneously use the first equation of motion (<i> = —V • B) to substitute for 
— $V • B in the action in order to replace it with $ 2 , since this would change the order of 
the action. 

It is interesting to examine the form of these dual theories. The potential for A/B 
essentially gets mapped into a non-canonical kinetic term for the F/$ (4-form/scalar) theory 
respectively. The canonical kinetic terms in the A/B picture, on the other hand, give rise 
to simple quadratic potential terms in the corresponding 4-form/scalar theories. This is 
important for several reasons. First, in this way an effective non-canonical scalar theory 
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arises from a very simple three-form theory. 4 Second, this immediately tells us that all 
scalar models dual to the three-form share the same simple quadratic potential. This is 
particularly important for standard slow- roll inflation 5 , where in the dual scalar picture the 
potential dominates over the kinetic terms. The fact that all models share the same dual 
scalar potential then implies that the form of the original three-form potential (which turns 
into a non-canonical kinetic term) is not important when computing, for example, the spectral 
index n s . We will return to this point later. 

Finally, one may wonder how one can come up with an effective single scalar theory 
dual to a three-form theory, which in principle possesses more physical degrees of freedom. 
However, starting with the most general canonical and minimally coupled 3-form action in 
4d, as we do here, guarantees that such a dual single scalar field description always exists. 
This is the case, because 1) the canonical kinetic term for the three-form dualises to a <3? 2 
potential 2) the three-form potential is a function of A 2 only, because in 4d this is the only 
covariant scalar combination that can be built from a 3-form - thus the potential only depends 
on one effective degree of freedom: A 2 and 3) these two degrees of freedom, <3? and A 2 , are 
related via an equation of motion, leaving only one effective independent degree of freedom, 
thus explaining the existence of a dual single scalar description. 

3.3 The perturbed action from scalar duality 

Now we wish to show that one can derive Eqs. (2.16) and (2.20) using the results above. In 
particular, we employ the expressions 

B 2 = -x 2 , X = -gTVpQVvQ = -12 2 B 2 (V{-6B 2 )) 2 , (3.11) 

which follow from the definition of the Hodge dual to A, Eq. (3.4), and from Eq. (3.6) 
respectively, to find 

x=v% md ^ = ^- (3 - 12) 

From the discussion of dualities, and in particular £q, we see that the P(X,4>) theory dual 
to our minimally coupled three-form theory was parametrised by 

P(X, 4>) = --cj) 2 - 12B 2 V'{-6B 2 ) - V(-6B 2 ) 



- 2 + xVx-V(Sx 2 ). (3.13) 



Differentiating this expression with respect to X and using (3.12) we find P t x, P,xx and 
P,xxx which upon substitution into 

c '= P ,? yP ' A = X 2 P t xx + \x z P^xxx , (3.14) 
P,X + 2A Pxx 3 

which are the expressions for c 2 and A for a P(X, <fi) theory [17], confirms the expression for 
the speed of sound (2.16) and for A (2.20) we found for the three-form theory. 

4 The reader might have noticed that there are differences in factors and signs when comparing (3.3) 
and (3.7) with corresponding equations presented in [24]. The version presented here corrects a small number 
of typographical mistakes in that work. 

5 By this we here mean e, r\... <C 1 as well as requiring no rapidly varying speed of sound. 
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This constitutes the result we are looking for. If we had only known about the perturbed 
action in terms of ( from a P(X, (f>) theory, we could have used the argument above to write 
that perturbed action solely in terms of three- form quantities, hence providing an alternative 
route to our earlier result. 

The power of this approach is that we can also readily probe higher order statistics, 
such as the trispectrum, since the quartic action has already been calculated for P(X, <p) 
theories [40, 41]. This involves the quantity II, defined in Ref. [41], which for completeness, 
we calculate to be 

tt _ 1 V ,X 'XXX J v ,x ,xxx 1 v ,x v >xxxx f"*1^ 

40 V? xx + 40 V xx 40 V xx ' ^ 0) 

in terms of the three-form quantity X- 

3.4 Confirmation of our gauge choice 

Finally we would like to confirm that the gauge we have chosen in §2.2 is equivalent to the 
gauge chosen when the perturbed action for the P(X, <j>) theory is calculated [17, 18]. There 
the choice is usually expressed as 5cp = 0. Utilising the expression 

X = 12 X 2 (V'(6 X 2 )) 2 , (3.16) 
and perturbing both sides, we find this leads to the expression 

a = 3xC-^ L ai, (3.17) 
v ,xx 

which is the same as the anti-symmetry condition (2.10) after we fixed the gauge choice 
a = a 2 xip- It is then possible to verify that the symmetry constraint on the three-form gives 
Eq. (2.8), and the other substitutions remain unchanged. In this way the equivalence of our 
gauge choice with the usual choice made in scalar field theories is confirmed, in contrast with 
the choice made in other works [27]. 



4 Non-Gaussianities 
4.1 Correlation functions 

We will now proceed to utilise our results thus far to compute observables for the three- 
form inflationary theory under consideration. Because our action is identical to that for 
a non-canonical scalar field (except with c 2 and A being expressed in terms of background 
three- form quantities), the calculation follows precisely that of Garriga & Mukhanov [12] to 
calculate the power spectrum, and Seery & Lidsey [17] and Chen et al. [18] to calculate 
the three-point function. Here we simply provide an overview of the main steps and the 
important results. 

The 2-point correlation function for the curvature perturbation £ is defined as 

(C(k 1 )C(k 2 ))=(2vr) 5 ( 5 3 (k 1 +k 2 )^, (4.1) 



and is calculated using the second order action Eq. (2.17). As found by [12] the end result is 



where * indicates that the expression is evaluated at horizon crossing c s k = aH. The spectral 
index n s is then given by 

l-n.=* + A + * (4 . 3) 



The slow-roll approximation (i.e. neglecting \ in the equation of motion for x) implies 



1 



that 

Xyv _ 
X 3 X 2 V" 2' 
Substituting this into the Friedmann equation we obtain [23] 



,X 2 I <■ 



H 

IP 



3 Vy 

-x— 
2 A V 



1 



Eliminating the term in brackets from the last two equations we can write 

X,N 4 1 V 2 



X 



9X 2 XV 



(4.4) 



(4.5) 



(4.6) 



■x 



This expression allows us to easily compute e from the definition of e to find e/eH = 2e+0(e 2 ). 
Similarly, we find that c s /c s H ~ 0(e 2 ) and also X/XH ~ 0(e 2 ). This means that to first 
order in the slow-roll parameters, the scalar spectral index is simply n s = 1 — 4e. 

The ratio of tensor to scalar perturbations was computed in Ref.[24] and shown to be 
related to the slow roll parameter e as 6 



16c.,e, 



(4.7) 



mirroring the analogous expression for non-canonical scalar field models. One can also com- 
pute the spectral index for tensor perturbations nt, finding rit = — 2e to first order in slow-roll 
parameters as usual. 

The third order action (2.19) is needed to calculate the three-point correlation function. 
First, however, it must be written in a suitable form, which follows from the use of the ip 
substitution, and is [17, 18] 



5 3 



dtd 6 x 



S I 



1 



+ 2A 



+ 



a 3 e 



3 + 3c 2 )CC 2 



ae , 
c{ 



2e s + l-c 2 s )adC,) 2 -2a-MdC,){da) 



+ 



a 3 e d 



rC + ^(dO(da)d'a + -(d\)(daY + 2/(0 



2a 



4a 



(4.8) 



where f{£)8L/8£\\ indicates terms proportional to the functional derivative of the Lagrangian 
evaluated at first order in Q. This would be zero if Q was Gaussian. Such terms can be removed 
by a field redefinition (though one must recall the redefinition when forming correlation 
functions so as to form corrections of Q itself). This form of the action thus identifies the 
relevant interaction vertices which contribute towards the three-point function. 

6 A more exact answer, fully taking into account that tensor and scalar modes freeze out at different times 
for models with c s 7^ 1, can be found in [42]. Here we just note that this extra effect means that the actual 
tensor-to-scalar ratio r is in fact smaller than naively expected from r = 16c s e for c s < 1. 



- 12 - 



At tree level in quantum field theory, and in the interaction picture, the In-In (equal 
time) three-point correlation function is given by the expression 

(C(t, ki)C(t, k 2 )C(t, k 3 )) = -i f dt'([C(t, kOC(i, k 2 )C(t, k 3 ), iMOl) , (4-9) 

where -ff; n t is the Hamiltonian evaluated at third order in the perturbations and follows 
directly from (4.8). Vacuum expectation values are evaluated with respect to the interacting 
vacuum By convention the 3-point correlation function is parametrised by the amplitude 
A. 

(C(k!)C(k 2 )C(k 3 )) = (2^) 7 5 3 (k! + k 2 + k^-^A , (4.10) 

where, again by convention the power spectrum P^ in the above formula is calculated for the 
mode K = k\ + /c 2 + /c 3 . Evaluating (4.9), one can determine A 

In principle A is a general function of the three Fourier modes, which are related by 
the condition ki + k 2 + k 3 = (within the slow-roll approximation the full form of A is 
given in Ref. [18]). For a given shape of non-Gaussianity, however, (see figure 3) the size 
of non-Gaussianity can be adequately characterised by a single- value measure /nl- For an 
equilateral shape (i.e. one peaking in the limit k\ ~ A; 2 ~ /c 3 ), this can be defined as [43] 

/ N f ' = 30 Ak ^ k f k3 , (4.11) 

where amplitudes are matched at ki = k<i = k% = K/3. Note that we follow the WMAP 
sign convention here, where positive /nl physically corresponds to negative-skewness for the 
temperature fluctuations. The parameters controlling the overall size of / N l Ui1 are c s ,S,A, 
which we now summarise for convenience: 

„2 _ V ,XXX v _ V ,X x_ V ,X V ,XXX (A 11\ 

Following [18] we may now compute in the slow-roll regime, finding the result at leading 
order to be 



4.2 Example I: A power-law potential 

We now consider some concrete examples, and begin by considering a three-form model 
self-interacting through a simple power-law potential 

C = -~F 2 -V A 2 p, (4.14) 

that is, V(A 2 ) = VqA 2p = Vq (6x 2 ) P , where p is a constant. 

This is a special example, where the equivalent P(X, 4>) theory is relatively simple, and 
is given by 

£ = (2p-l) - {—A (4.15) 



VqJ \2Ap 2 J 2 



7 Note that (4.13) differs from the expression found by Chen et al. [18] by an overall sign. This is because 
Chen et al. use a sign convention opposite to that used by WMAP and throughout this paper. 
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where one now sees explicitly that the 3-form potential has been mapped into a non-canonical 
kinetic term for the effective scalar (f>. As long as 2p — 1 7^ 0, one can check that equa- 
tions (2.16), (2.18) and (3.14) equivalently yield 

4 = 2p - 1 , 

-= 11 ^, (4.16) 

5 3 2p - 1 ' K J 

and we can immediately make observational predictions for the theory, finding several inter- 
esting results. 

Here both the speed of sound and A/E are constant, and we have n s — 1 = — 4e. 
The constancy of the speed of sound has the interesting consequence that the observational 
requirement of obtaining a (near)-scale-invariant spectral index n s forces e to be close to zero 
and slowly- varying here (compare (4.3) and also [43]). For this power-law model, slow-roll 
therefore becomes an observational requirement in contrast to models with varying speed of 
sound. 

We now need to calculate this quantity N e-folds before the end of inflation. It was 
shown in [23] that for a power law potential, the value of the field at this time can be 
estimated to be 

2 _ 2 4 _J_ (A 17^ 

Xn ~3~ 18plT2iV' [ ' 

which upon substitution in (4.5) gives 

W «> (4.18) 

Assuming that N ~ 60 is required to solve the horizon problem, we predict that the spectral 
index on observationally relevant scales for a three-form with power law potential is 

n s « 0.97, (4.19) 

and independent of the value of the power p. Naturally the spectral index will be closer to 
scale-invariant the longer inflation lasts. Substituting for and A/E into (4.13) we obtain 
the dependence of on <? s illustrated in Fig. 1. It clearly shows that a small speed of 

sound leads to a large Z^" 1 as expected. Substituting for the speed of sound for a power law 
potential in Eqs.(4.7) and (4.13) we can relate and r. Figure 2 illustrates how these two 
quantities are related, and shows the region of the parameter space (r, /jvx) that is allowed 
given current bounds. 

Phenomenologically we find that no large non-Gaussian amplitude of the enfolded or 
orthogonal types (which have peaks in the folded limit k\ ~ k<i ~ 2/C3) can be generated here 
in contrast to generic single field inflation models. To see why, it is useful to notice that one 
may express A/S as [17, 43] 

> ifvx+i ;\ (4 20) 



E 6 V c? 



where 

/X = 3^' £s = i^? eX = ~lPdX- (421) 
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10 4 0.001 0.01 0.1 1 5 

Figure 1. Dependence of on c 2 s for the power law potential V = Vo(6x 2 ) p and N = 60. A 

large and generically negative non-Gaussian amplitude is found. The shaded region is disallowed by 
the WMAP 2a bound -214 < f^ il < 266.[15] We recall that c 2 s = 2p - 1 for this model. 

r 




Figure 2. The solid lines show how the parameters and r are related to each other for the 

power law potential, V = Vo(6x 2 ) p and for N = 50, 60, 80 from top to bottom. The shaded region is 
disallowed by the WMAP 2a bound -214 < f^ il < 266. [15] 



Now, since the speed of sound is constant for the power-law model, e s , fx are identically zero 
here. A large and predominantly orthogonal (or enfolded) amplitude, however, requires a neg- 
ative and non-zero fx (assuming positive e and < c s < 1)[44] 8 . As such the non-Gaussian 

8 In general P(X, inflation models a predominantly orthogonal or folded shape can be generated by finely 
balancing the contributions from the <^" 3 interaction vertex (which depends on A/E) against the other vertices 
such that the generically predominant equilateral shape contributions cancel out. In this way the otherwise 
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Figure 3. The dimensionless bispectrum A(k\, fe, ksj/kikzka for the power law potential, V = 
Vb(6x 2 ) p with p = 0.505, corresponding to c s = 0.1, and with N = 60. Triangular shapes denote the 
equilateral, squeezed/local and enfolded limit clockwise from top right. A predominantly equilateral 
shape is found. 

shape found here is always predominantly equilateral. Furthermore, since parameters A and 
S are simply related to c s via (4.20), the non-Gaussian amplitude is completely controlled 
by (? s for a power- law potential. This means that /^i!" 1 ~ ^( c 7 2 ) f° r subluminal and posi- 
tive speed of sound c s , i.e. p e [1/2, 1]. Figure 3 illustrates these points. Finally note that, 
analogous to e.g. DBI inflation, any sizable level of non-Gaussianity has a negative f^ 1 ^ 
associated with it here. 

4.3 Example II: Exponential potential 

Let us now consider an alternative potential. Even though it is fairly simple in its analytical 
form, the exponential potential, V = Vbexp(/3A 2 ) = Vq e~x.p(6f3x 2 ), becomes very complex 
when written P(X, (j>) theory, 



where W(x) is the Lambert-W function and x = X/12(3Vq. Dealing with this model in 
the P(X,(ft) description is therefore very difficult. This is a good example showing that 
performing the calculations in the original three-form theory is far simpler than going to the 
P(X, (f)) description. 

Inflation ends when e ~ 1 which for this potential takes place when 



subdominant orthogonal or folded configurations are brought out. As stated above, this delicate cancellation 
relies on /x-dependent contributions though, whereas f x = identically for the power-law example given. 




(4.22) 




(4.23) 
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This expression implies that inflation only ends if > 1/3. In particular, if (5 S> 1/3 we can 
approximate 

In this case, the solution of (4.4) and (4.5) is non-analytical and we must resort to an 
approximation. Defining x 2 = 2/3 — y, since we know that inflation occurs very close to 2/3, 
we obtain 

xl = I ~ W TTTSv' (4 ' 25) 

The slow-roll parameter N e-folds before the end of inflation is then e^v ~ 1/(1 + VQN). Here, 
c 2 = 1 + 12/3x 2 and it can be verified that c s /c s H = -e 2 /x 2 ~ 0(e 2 ), hence, for AT = 60, we 
obtain for this model 

n s « 0.97, (4.26) 

independent of the value of /3. The fact that the value of the scalar spectral index is indepen- 
dent of parameter j3, or p in the case of the power law potential, and we obtain the same value 
in both examples, should not come as a surprise. The reason becomes clear when considering 
the dual P(X, (f>) theory, which for any three-form potential, has the same quadratic scalar 
potential. Since we are in a slow-roll regime, the functional form of the kinetic term is not 
important with respect to the potential which leads to identical results for the spectral index. 

It was found above that the choice of the exponential three-form potential offers an exit 
from inflation only if /3 > 1/3. Taking x% ~ 2/3, this puts the lower bound, c 2 s > 11/3, 
which means that the speed of sound is superluminous. 

The amplitude of the three-point function is controlled by c 2 and A/X, which are given 

=1+1**-, * jgva+jftfl . (4 . 27) 

Varying j3 and substituting for c 2 and A/S N e-folds before the end of inflation into (4.13), 
we obtain the dependence of f^ tl on c 2 which we show in Fig. 4 for N = 60. In Fig. 5 we 
see how r and /nj™' relate to each other. In particular, we observe that the large values of c 2 
render this model disfavoured by current bounds on the ratio of tensor to scalar perturbations, 
r < 0.2 for N = 60. 9 Of course, if we allow for more e-folds of inflation this reduces the 
relevant e and hence r. A crude estimate for the minimal amount of inflation to bring this 
model into agreement with current la bounds on r may be obtained by assuming a speed 
of sound right at the lower bound of c s > ^11/3 and r ~ 0.2, yielding N > 62 for the 
exponential model considered here. Finally we see in Fig. 6 that the sign of /mt" ~ 0.2 is 
positive, albeit with a rather small amplitude. The non-Gaussian shape here is predominantly 
equilateral, but also picks up contributions in the enfolded/orthogonal limit. 10 



5 Conclusions 



In this paper, we have considered three-form theories of inflation and their observational 
footprints. The salient results may be summarised as follows. 

9 Note that this bound significantly depends on the assumption of no running spectral index n s . If a running 
n s is allowed the bound weakens to r < 0.49. 

10 Note that considering departure from slow-roll could modify this statement, e.g. potentially rendering the 
shape predominantly enfolded or orthogonal [45, 46]. For further details on fast-roll corrections also see[47, 48]. 
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Figure 4. Dependence of f^™ 1 on c 2 . for the exponential potential, V = Vbexp(6/3x 2 ) and N = 60. 
A small and generically positive non-Gaussian amplitude is found. Notice the lower bound c 2 . > 11/3 
for this model. 

r 




Figure 5. The solid lines show how the parameters and r are related to each other for 

the exponential potential, V = Vbexp(6/?x 2 ) and for N — 50,60,80 from top to bottom. The 
shaded region is disallowed by the WMAP 2a bound r < 0.2. [15]. End points for the solid lines at 
ff 1 ~ 0.215 correspond to the lower bound c 2 s > 11/3 for this model. 




• We have calculated the action for a three-form theory perturbed about a FRW back- 
ground, isolating the perturbed quadratic and cubic action expressed purely in terms 
of the comoving curvature perturbation £, as well as background and three- form quan- 
tities. 

• We then reviewed dualities between various p-form models in four dimensional space- 
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Figure 6. The dimensionless bispectrum A(k%, &2, k^/kxkzkz in the slow-roll limit for the expo- 
nential potential, V = Vo exp(6/?x 2 ) with (i — 1, corresponding to c s ~ 3 and for N = 60. Triangular 
shapes denote the equilateral, squeezed/local and enfolded limit clockwise from top right, ~ 0.2 

and a predominantly equilateral shape is found. 

times. In particular, these allow the linking of the three-form theory under considera- 
tion to an effective scalar description with a P(X, 4>) action, due to the fact that the 
most general canonical and minimally coupled 3-form action in 4d has only one indepen- 
dent degree of freedom. Since the three-form potential is mapped to the scalar field's 
kinetic term in a non-trivial way, such setups generically result in highly non-canonical 
P(X, 4>) theories, which do, however, all share a simple <E> 2 potential. 

• Despite the complexity of the mapping, we were, however, able to use this duality as 
a formal tool to give an alternative derivation for the perturbed quadratic and cubic 
action. To do so we used existing "non-Gaussian machinery" for P(X, 4>) theories, but 
translated the end-result back to three- form quantities (essentially re-expressing c s , £ 
and A in terms of the three- form variable x)- We showed that the result is identical to 
that obtained via the direct three-form computation. 

• Finally, using the tools developed, we explored the observational and particularly non- 
Gaussian features of two example three- form models: 

— First, we investigated three-form inflation with a power-law potential. This simple 
setup produces a constant speed of sound and generically equilateral non-Gaussian 
shapes with negative f*rrf\ similar to DBI-type inflation. Observational limits on 
this form of non-Gaussianity constrain this model. 

— Secondly, we explored the phenomenology for a three-form inflationary model 
with an exponential potential. We found the non-Gaussianity was once again 
predominantly of equilateral shape, but in this case was unobservably small. We 
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found, however, that the model can be constrained from limits on the tensor-to- 
scalar ratio r. 

• An interesting general result is that the spectral index for all three-form inflationary 
models is n s ~ 0.97 to leading order in slow-roll when 60 e-folds of inflation occur. This 
was shown explicitly for the two models at hand, but will hold much more generally. 
The value of n s is therefore uniquely predicted once N, the number of e-folds of inflation, 
is specified, independent of the exact form of the three-form potential and as long as 
e < 1. 

In conclusion, three-form inflation has a number of interesting features. Simple three- 
form potentials lead to interesting models of inflation with potentially large non-Gaussian 
signatures, but with a spectral index which is within the preferred WMAP bounds. It 
would be interesting to extend our studies of explicit three-form potentials to more complex 
potentials which have no analytic form for the dual scalar theory; to use our framework to 
compute the trispectrum; to consider a multi three-form action. We defer these questions to 
future studies. 
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A Appendix I: Dual fields 

Here we give some pedagogical detail on Hodge-dualising p-forms in n dimensions. We begin 
by defining a totally anti-symmetric tensor e as 

j ( +1 if (a±, Q2-..a p ) is an even permutation of (1, ...d) 

e ai,a 2 ...a p _ < — 1 if Q 2 ...a p ) is an odd permutation of (1, ...d) (A.l) 

V Isl [ otherwise, 

where d is the number of spatial dimensions and we define dt to be the number of temporal 
dimensions, so d = 3 and dt = 1 for the cases we will consider in this paper which have 
signature ( — h ++). Note that for a diagonal metric consequently (— l) dt = sign(g), where 
g is the determinant of g^ u . Lowering all indices with the metric g^ v (note this is a valid 
procedure since we are explicitly dealing with the tensor e, not its associated tensor density) 
one finds 

f +1 if (ai, a.2---Otp) is an even permutation of (1, ...d) 
—1 if (ai, ct2---a p ) is an odd permutation of (1, ...d) (A. 2) 
otherwise. 

A particularly useful identity we use repeatedly is 

e ai .... an e^ = {-) d H\8f^ 2 ...8^. (A.3) 



e ai ,a 2 ...a p = (-) dt VW\< 
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Let us now consider an arbitrary p-form living in an n-dimensional space, where we choose 
a coordinate basis and write 

P = -<P ai ,... a Jx ai A ...dx a ", (A.4) 
pi 

In terms of the totally antisymmetric tensor e the Hodge (*) dual of this p-form is given by 

We recall this means any p-form has a dual which is a (d — p)-form. In particular the 
three-form A and four- form F that make up (3.3) can therefore be expressed as (3.4). 
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